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ABSTRACT 

This thesis presents a microcomputer based, 
computer-aided design program -for tactical military 
tropospheric scatter radio systems. The program has the 
capability of predicting the system performance and 
reliability for both analog (FM/FDIi) and digital 
trap osc at ter radiolinks. Propagation gain generated bv 
elevated tropospheric ducting is called height gain. A 
height gain computational model for specific elevated 
tropospheric ducts is derived from statistical radiosonde 
data. A terrain profile plot, real-time radiosonde data 
analysis, and the probability of error for digital 



rad iolinks 



are 



provided as program options. 
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I. INTRODUCTION 



A. PURPOSE 

The development of digital transmission techniques in 
military commun i cat i ons systems has generated a renewed 
interest in tropospheric scatter- Digital troposcatter 
radio systems have required changes in the present 
planning and engineering design methods. The military 
design engineer has the responsi b i 1 i t y to predict system 
performance and reliability within limited planning 

time. Under certain tactical situations, the urgency to 
restore vital commun i cat i ons has forced the engineer to 
provide several alternative designs for a single 
communications link- These requirements have justified 
the need for an efficient, computer-aided design method 
to accomplish the task. The purpose of this thesis 
research is to develop a tropospheric scatter 

communications system design program. The objective is 
to design a mi crocomputer-based program which will 
predict the path reliability and system performance for 
troposcatter radiolinks. This program will introduce 

1 • : 



the 



manual 1 y generated , 



several improvements over 
nomograph-based design methods 
program will estimate the height 
system antenna height -for specific 
ducts. Other features include a 
radiosonde weather data analysis, 
the probabilty of error for binary 
signaling techni ques . 



Speci f i cal 1 y , the 
gain at the average 
el evated tropospher i c 
terrain profile plot, 
and the calculation of 
and orthogonal digital 



B. BACKGROUND 

Tropospheric scatter communications systems propagate 
microwave energy beyond 1 i ne-of -si gh t (LOS) or “over the 
horizon" by taking advantage of the refraction and 
reflection phenomena in a section of the earth's 
atmosphere called the troposphere. Typical military 
troposcatter systems use transmitter power outputs from 1 
to 10 kW, parabolic type antennas, and sensitive 
broadband Fli receivers with front-end noise figures (NF) 
between 2.0 - 4.0 dB. These systems are most often 
employed in military applications because of the limited 
available channel capacity [Ref. l:pp. 253— 254 D. Several 
system advantages are summarized as follows: 

1. Fewer terminal and/or relay stations are required 
to cover the same transmission path as compared 
to tactical microwave 1 i ne-of -si ght radiolinks. 



2. Reliable multichannel communications can be 
installed across long distances of hostile or 
inaccessible terrain. 

3. Standard transmission ranges are suitable for the 
tactical military field environmemt for radiolinks 
from 30 to 200 miles. 

A fundamental design parameter for troposcatter links 
is the allocation of allowable channel degradation to 
each link making up the total communications system. In 
analog scatter systems, the quality parameters are system 
availability and signal-to-noise ratio in the voice 
channel. However with the deployment of digital 

troposcatter radio systems the performance of the 

digitized voice channels under fading conditions is 
si gni f i can t 1 y different than the performance of analog 
voice systems under the same conditions. The primary 
measure of transmission quality for a digital system is 
its error performance. There are two main sources of 
transmission errors: (a) long-term error rate which will 

occur because of equipment degradation, channel 

interference, and long-term power fading; and (b) 

multipath fading [Ref. 2s pp. 4-5 D. 



c. 



RELATED WORK 



1 . Performance Models for Troposcatter Links 

Monsen CRef. 33 has derived performance models 
for analog FDM/FM and digital quadrature phase shift 
keying (QPSK) systems on troposcatter communications 
links. The analysis included the effects of signal level 
variations and multipath delay dispersion. His 

performance criterion was the outage probability rather 
than either average error rate or median si gnal -t o-noi se 
ratio (SNR) values. Outage probability performance in a 
digital system was derived for two different modems: the 

Decision-Feedback Equalizer (DFE) and a transmitter 
time-gating technique, the Distortion Adaptive Receiver 
(DAR) . Analog system performance was determined for the 
percentage ' of the time that the voice channel SNR 
including thermal noise and multipath delay dispersion 
effects were below a specified threshold. 

The performance test results for a Distortion 
Adaptive Receiver (DAR) was presented in a paper by 
Zawislan [Ref. 43. The critical problem with digital 
troposcatter is intersymbol distortion produced by the 
multipath dispersion of the fading channel. To resolve 

the distortion adaptive receiver 

J. 3 



this problem, 



was 



implemented. The DAR modem employed QF'SK modulation with 
adaptive matched filter demodulation. A complete 

functional description was explained in the reference. 
The advantage to this approach was that it did not 
require equalization at the receiver to correct 

intersymbol interference and it provided near optimum 
performance using the adaptive matched filter receiver. 

Typical bit error rate performance for the DAR 
on a troposcatter channel is shown in Figure 1-1 for 
different values of rms multipath delay dispersion. For 
low multipath spread , the fading on each diversity 
channel has a Rayleigh distribution and the per f crmance 
is indicated by the "flat fading" marked curve. Note 
that the performance improves rapidly with increasing 
multipath spread and, at a 10 error rate a 10 dB 
improvement is achieved for multipath curve CBU. However, 
as the multipath spread increases (with increased path 
length), a irreducible error rate occurs as shown in 
curve LC1 due to the transmitter time-gate limitation. 
To compensate for this problem, the DAR is modified by- 
utilizing a dual frequency pulse waveform. This 

technique is also used in the modem employed by the 
AN/TRC-170 tactical digital troposcatter radio set. 



SINGLE PULSE DAR PERFORMANCE 




EB/NO (P6) 



(DUAL EXPLICIT DIVERSITY) 



Figure 1—1 Distortion Adaptive Receiver Performance 
(After Zawislan, Ref. 4:p. 3) 
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2. Tropospheric Elevated Duct Models 



Several mathematical models have been developed 
to calculate the ducted signal -field strength based on 
the re-fracti vi ty profile and by modeling the tropospheric 
duct as a laterally homogeneous-layer atmospheric 
waveguide [Ref. 5 : pp - 22] . A detailed mathematical 

analysis was presented by Marcus and Stuart [Ref. 61 . 
These authors have developed an effective modal solution 
computer program, called "DUCT", which requires a large 
processing capabilty of a main-frame computer [Ref. 71 . A 
version of this program, called "PDUCT n , is available at 
the Naval Postgraduate School Computer Center [Ref. S] 
and [Ref. 5:pp. 23-24]. 

Weston [Ref. 5:pp. 58-70] has proposed a 

statistical model, based on "PDUCT" predictions for 
selected ducts of interest. The height gain data, for 
these ducts, was used to derive model coefficients. 
These unique coefficients were stored according to their 
respective transmitter and receiver heights in a matrix 
form. Thus, given a receiver height and range from the 
transmitter for a selected duct /frequency data point, the 
coefficients for each receiver height are obtained and 
used to produce the correspond! ng height gain curve. 



This approach required a significant number of main-frame 
computer "PDUCT" computations to establish height gain 
data for all the ducts of interest. 

In December 19S3, a mi crocomputer-based pro- 
gram, called " MINI DUCT “ Version 1-1, based on Knorr's 
mathematical model was developed by Nagel, CRef. 91, at 
the Naval Postgraduate School to calculate elevated 
ducted signal levels at various frequencies- This 

program used either historical radiosonde data or current 
elevated duct information CRef. 101- The radiosonde data 
represented the atmospheric temperature, pressure, and 
humidity at specific altitudes over a five year recording 
period- This model was valid for the case where the 
transmitter is at the optimum duct coupling height and 
the receiver is located either below, within or above the 
duct- The optimum coupling height is the altitude above 
the surface where the gradient of the modified refractive 
index becomes negitive. 



II. TR0P05HERIC PERFORMANCE CONS I DERAT I QN5 



A- SCATTER PROPAGATION 
1 - General 

At frequencies above 30 MHz three propagation 
mechanisms can carry energy beyond the horizon: (a) 

variations in the refractive index in the troposphere can 
scatter radio energy, <b) horizontal 1 y-strat i f i ed abrupt 
changes in the refractive index can cause reflection, ana 
(c) atmospheric regions of negitive modified refractive 
index gradients can introduce ducting. Refractive index 
and tropospheric ducting will be thoroughly discussed in 
the following section. The forward scattering of radio 
signals is the most dominant propagation mechanism at the 
frequencies of 0.3 to 10 GHz. [Ref- ll:p« II 
2. Received Scattered Field 

The index of refraction depends on pressure, 
humibitv, and temperature. Slight variations in these 
quanitities, caused by atmospheric turbulence, will 



produce random f lucuations in the refractive index. When 
" ' ' ' ’ 1 *“ thi s 



an 



el ec tromagnet :i c 



wave 



propagates 



through 



inhomogeneous medium, energy will be scattered out from 

the original incident direction. The 

turbulent-scattering theory, CRef. 12:p. 3453, has shown 

to a first approximation that the index of refraction 
fluctuations can be replaced by a model of so-called 
“blobs", inhomogeneities of different dielectric 
constants randomly distributed. If these “blobs" are in 
the common volume formed by the transmitter and receiver 
antenna beams, the complex received field can be 

described by CRef. 13:p. 1461: 

Re^- f>. 

1 i r>- i ^ 



where* m is the 

vol ume , A is 
J 



nuimber of “blobs" in the scattering 

the amplitude, and (p is the phase of a 

J 



wave scattered by a single "blob". Assume m to be very 
large and the blobs are spherical and uniformly 



distributed thr ouigh 


the 


scatter i ng 


vol ume . 


Then 


the 


phase difference 


of 


the waves 


scatter ed 


by 


the 



i nhomogenei t l es at the top and bottom of the volume will 
be CRef. 13:p. 1471: 



A4> 



4 n h sin y 
A 



( 2 - 2 ) 
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where y = transmitter, receiver antenna elevation angle 



h = mean thickness of the scatter volume (m) 

X = transmitter wavelength (m) 

The variables in Equation 2-2 are described in Figure 2-i. 
By observing the volume geometry, it can be determined 
that more "blobs" exist in the central volume region than 
at the top or bottom. This indicates a non-uniform 
"blob" distribution. 

Many phase variations wi IT occur over many 
phase cycles of length 2n . The antenna elevation angle 
will be less than 2 degrees for most systems and the 
ratio of scatter volume thickness to operating wavelength 
will be very large. These conditions make the scattered 
phases uniformly distributed from zero to 2n 

[Ref. 1 3 : p . 148], 

The amplitude distribution will now be deter- 
mined. The real and imaginary components of the complex 
received field, Equation 2—1, were resolved into two 
random variables X and V. The number of blobs, m, is 
assumed to be large, so by the Central Limit Theorem, 
both X and V will approach a normal "gaussian 11 
distribution. Finally it was shown that these variables 




Figure 2—1 Scatter Volume Geometry 
(After Beckmann, Ref. ll:p. 147) 
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are also uncorrelated and independent. 



Then 



thei r 



two-dimensional normal distribution can be transformed 
into the Rayleigh distribution. 

Turbulence, mixing, and wind continuously 
change the positions and structure of the inhomogeneities 
making the random terms of Equation 2-1 functions of 
time. The reality of nature is the inhomogeneities do 
not vary isotropically because the index of refraction 
changes more rapidly with altitude than with distance. 
There can also be a str at i f i cat i on or layering of 
inhomogeneities near a particular height (or heights) 
that changes the assumed spherical shape of the so-called 
"blobs". Under these conditions the uniform phase 
distribution may change to an unknown di stri but i on. 



B. TROPOSPHERIC 
The index of 
the ratio of the 
(EM) radiation 
con veni ent parameter 
[Ref . 1 1 : p. 14D : 



n , of air is 
in a vacuum of 
vel oci ty in 
is ref r act i vi ty , 



DUCTING, 
r ef r act i on , 
vel oci t y 
to the 



defined as 
el ectromagnet 
the medium. 

N, defined 



1 c 
A 
as 



N = ( n - 1 ) X 10 6 



(2-3) 






Another parameter called modified refractivi ty, 



1 s 



defined as CRef. 10:p.9D: 



Mi 



M = N + 0 . 15 7 h 



(2-4) 



where h is the altitude in meters above the surface. The 
modified refractivity accounts for the curvature of the 
earth so the presence of ducting can be easily determined 
by observing the M-gradient on the M versus height plot. 

Refraction of incident radio waves a cros s a 
discontinuity of refractivity is described by the 
principles of Snell's Law. It is important to remember 
that the wave "bends’ 1 towards the higher value of 
refracti vity , and the more dense a material the higher 
its n. Since the density of the atmosphere decreases 
with height, we expect that a wave will bend back 
downward from a geometric straight path. Whenever the 
refractive index decreases sharply with height, 
radiowaves can be trapped and experience low-loss 
propagation for long distances. This condition is known 
as tropospheric ducting CRef. 1 1 : p . 291. 

The following conditions must be satisfied for a duct 
to occur: (1) the modified refractive index gradient 
shall be equal to or more negitive than 0 M— units/km, and 



(2) this gradient should continue over a height of many 

wavelengths. The important duct parameters are the duct 

thickness, D, the intensity, M, and the optimum coupling 

height, H , A piecewise linear appr ok i mat i on to the 
opt 

modified refractivity (tr i -1 i near ) profile for several 
types of ducts are shown in Figure 2—2- There are three 
types of ducts: (1) surface or ground-based ducts, also 

called evaporation ducts when formed over water. Figure 
2-2a; (2) surface-based ducts from elevated refractive 

layers. Figure 2-2b ; and (3) elevated ducts from elevated 
refractive layers, Figure (s) 2 -2c and d. Note that all 

positive M-gradients are assumed at 118 M— units/km which 
corresponds to a standard atmosphere. Once the slopes 
are identified, the important duct parameters are quickly 
determined [Ref * 9:pp, S— 123. 

Tropospheric ducts more often occur as ground-based 
ducts because of both. evaporation and advecticn. 
Evaporation of water vapor from the surface of the sea 
may create a zone of high humidity (i.e. high refractive 
index; below a r e g i on of d r i e r air, Advection , defined 
as the movement of one air type over another, may cause 
hot dry air (from the land) to be blown over cold wet 
air, producing a region of low refractive index above a 
region of high refractive index- Such a duct may also 



Elevot ion , 




(o) Ground-based layer, (b) Elevaled loyer, 
Ground-bosed aucl Ground-based due! 

a 




( c ) E levoied loyer, 
Elevoied duel 



( d ) E levof ed loyer and 

dud approximated with 
on M - discontinue y . 



Figure 2-2 Tropospheric Duct Descriptions 
(After Marcus, Ref. 7:p. 2-5) 




form when warm dry air is blown over cold ground. 
Radiation cooling can produce temperature gradients which 
cause ground-based ducts- Air next to the ground becomes 
cooler and the duct becomes thicker as the night 
conti nues- 

When morning solar heating warms the air next to the 
ground, a region of rapid decrease of refractive index 
with height produces an elevated duct- However these 
ducts quickly disappear because continued ground heating 
increases convection mixing and destroys the stable 
elevated layer- Elevated ducts may form for several days 
by a subsidence inversion- Hot air rises at the center 
of a high pressure region and spreads out horizontally, 
cooling as it slowly descends. This produces a boundary 
with the slightly colder air near the surface. The 
increasing temperature with height at the boundary forms 
the subsidence ( temperature) inversion- Changes in 
temperature and/or humidity may cause related changes in 
the refractive index within the boundary interval 
CEef- llspp. 33—35 1 - 



Figures 2— 3a thru d 
within ground— based ducts, 
at elevation angles close 
surface , 



i 1 1 List rate propagat i ng ravs 
Rays leaving the transmitter 
to the horizontal will parallel 
departing ravs 



the earth ' s 



while other 



wi 11 




Figure 2-3 Propagation Effects of Ducts 
(After Hall, Ref. 9:p. 31) 
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either travel upward or downward. It the M-gradient is 
negitive and assumed uniform within the duct, rays 
beginning at elevation angles below the critical angle, 
defined as CRef. lisp. 31D: 



- >< 



(2 Am x 



10 ~ 6 ) 



(2-5) 



will strike the surface, and those above the critical 

angle will leave the duct. If ground reflection is 

neglected, Figure 2— 3a , a shadow region depicts minimum 

propagated energy. In the case where surface 

reflectivity is high, Figure 2-3b , the rays launched 

within the critical angle range will travel beyond the 

radio horizon. Figure 2-3c decribes the situation when 

the M-gradient is not constant within the duct. If a 

transmitter is positioned in the height region h , then 

a 



- 6 , 



rays within the + ^ (2 A M ^ X 10^) range will 
remain in the duct with bounces as depicted in Figure 
2-3b. Similarly rays transmitted within 



t he 1 0 2 = “ 4 (2 A M 2 X 10 range are trapped within the 
height range h . Finally a sur f ace— b ased duct from an 



b 

elevated 1 ayer 
Figure 2— 3d . In 
place at the top 



over rough terrain is 
this case the refractive 
of the duct [Ref. llspp. 



illustrated in 
“bendi ng '* takes 
30-322 . 
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For elevated ducts -formed by advection or subsidence, 

the position of the transmitter and receiver relative to 

the optimum coupling height, H ^will influence the 

opt 

propagation effects- The means in which energy enters or 
leaves an elevated duct can be described by the duct 
acting as a "leaky" waveguide,. Energy is "leaked" or 
coupled into the duct from the transmitter, and "leaks" 
out as the energy propagates along the duct 
[Ref- lisp. 36 1 . 

Because of the non-permanent char ac t er i st i c of 
elevated ducts, their effects are seldom an influence to 
troposcatter links, especially if they form above the 
common scatter volume. But the presence of tropospheric 
ducts can degrade the overall performance of troposcatter 
systems by changing the predicted transmission loss- The 
term that will change the total path loss is the duct's 
height gain, which is derived in Chapter IV. 



C . MULT I PATH CONS I DERAT I GNS 

The multipath fading model for a tropospheric scatter 
channel produces received signal fading- The received 
signal consists of the sum of a large number of 
time-variant, complex vectors having amplitudes and 
phases. The fading is caused by r an damn 1 y time-variant 



phases variations. At times the received signal vectors 
add destructively to decrease the mean received signal 
amplitude. While at other times, the vectors add 
construct i vel y , so the received signal is large. Thus 
the amplitude variations, or signal fading, are due to 
the multipath character i st i cs of the tropospheric 
channel . The channel can be modeled as a zero mean, 
complex-valued gaussian process, with the envelope of the 
instantaneous signal level being Rayleigh-distributed. 
This Rayleigh fading channel describes the short-term 
fading. When there are fixed constant regions of 
refractivity or stratified refractive layers in the 
vicinity of the common scatter volume the channel cannot 
be modeled as having a zero mean. In this case, the 
Rayleigh distribution does not apply and the channel 
approaches a Rice distribution [Ref. I4:pp. 456—458 j . All 
performance predictions in this study will assume 
Rayleigh 11 short-term " fading as the channel model. 

Channel performance will be degraded during periods 
of severe multipath fading. With digital systems, the 
voice user is unaware of any increase in background noise 
until the PCM (Pulse Code Modulation) outage threshold is 
broken. Once this threshold is passed the complete 
multichannel circuit will be unusable due to noise. The 



character i st i cs of fade outages for a typical digital 

troposcatter circuit can experience three (3) primary 

categories of fade outage. The first category occurs 

when the voice user is subjected to a single fade outage 

of duration less than 200 mi 1 1 i seconds. This outage will 

be hardly noticed. The second category of outage will 

have an outage duration ranging from 1/5 second to 5 

seconds. The user will detect this distortion but will 

continue to communicate following the outage. When a 

recurrence of short duration outages take place (e.g. 2 

to 4 outages per minute) annoyance rather than total 

disruption will occur. The third category are fade 

outages that exceed a subjective level of user patience. 

The Defense Communications Agency (DCA) [Ref. 15 J has 

specified five (5) ranges of fade outage conditions, 

refer to Table I . DCA has defined the fade outage in 

terms of a diversity signal-to-noise ratio threshold 

-4 

correspondi ng to a 10 bit error probability. 

Techniques used to counter multipath propagation are 
frequency diversity, space diversity, amplitude 

equal i z at i on , and channel equal i z at i on . If a modulated 
signal is simultaneously transmitted over the same 
troposcatter radiolink on two or more frequencies, the 
correlation between the individual received signals will 
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TABLE I 



VOICE PERFORMANCE CHARACTERISTICS 

FOR 

DCS TROPOSCATTER LINKS 



Outage Range Criteria 



Outage Probabi 1 i ty 



I 


See 


Note 


See Note 
-4 


1 1 


0.2 


sec.< Outage <5 sec. 


7.5X10-5 


III 


5 


sec.< Outage <1 min. 


7.5X10-3 


IV 


o 


< Outages/min. < 5 


2.5X10-4 


V 


5 


< Outages/min. 


1.0X10 



NOTE: 

Range Voice Performance Description 



I 

I I 

I I I 

IV 

V 



Outages with adequate fade margin. 
Outages with adequate fade margin 
and high frequency of occurrence. 
Call disruption possible. 

Marginal fade margin. 

Unavai 1 ab i 1 i t y of circuit. 



be small- This method o-F signal diversity is called 
frequency diversity. The important advantage of 
frequency diversity is that it requires a single antenna 
at each site- But the need for additional frequencies 
can increase the probability of co-channel interference 
among other operating transmitters. 

Uncorrelated short-term fading can also be achieved 
by separating the receiving antennas in space. This is 
known as space diversity- Horizontal and vertical 
polarization can ‘be used to distinguish between two 
space-separated signals- However horizontal and vertical 
polarization do not provide a satisfactory degree of 
noncorrelation of signal fading for efficient diversity. 

The multiplicity of signals provided by these 
diversity methods must be combined. The 
di versi ty— combi ni ng techniques are classified into: (1) 
selection or switching; (2) combining a desirable 
weighted combination of received available signals; and 
(3) a combination of selection and combining. 

In the selection techniques the di versi tv channels 
are scanned until one is found whose level exceeds a 
selected threshold- This may not necessarily select the 
best available signal. In the combining techniques, all 



diversity channels are simultaneously monitored and 
equally weighted- This is called equal-gain combining. 
In max i mal -rat i o combining, the weighting -factor o-f each 
channel is automatically adjusted to yeild the maximum 
si gnal -to-noi se ratio -for the total o-f all the diversity 
channels CRe-f. 12: pp. 453-454 U . 

Amplitude equalizers are designed to properly 
equalize the propagation channel -for minimum phase 
■fading. Channel equalizers balance the channel -for 
amplitude and mutlipath delay distrotion. They are 
typically adaptive transversal equalizers that consist ot 
tapped delay lines (TDL) with tap-weight multipliers 
(i.e. , amplifiers or attenuators); and control circuitry 
that adaptively vary the tap-weights in response to 



temporal channel variations CRe-f. 16:p. 



1 1 — 121 ! . 



III. TRQPQSCATTER SYSTEM DESIGN PROGRAM 



A. GENERAL PROGRAM ALGORITHM 

The main program, named "TROPQ" , i s presented as 
Appendix B. Figure 3-1 illustrates the program flow and 
the primary computational program modules. The program 
development and mathematical approach for the Radiosonde 
Data Analysis and Height Gain modules are described in 
Chapter IV. The remaining modules are formulated in this 
chapter. A program tutorial and compilation instructions 
are contained within the program. 

B. PREDICTION OF PATH LOSS 
1 . Gener al 

The basic median transmission loss will be the 
sum of several additive losses, expressed in decibels, 
[Ref. 171 z 



L 



T 



L + L , + L + L + L — G 
s d c a w t 



G + HG (3-1) 
r 



where L = free-space propagati on/scatter loss (dB) 
s 

L = knife-edge diffraction loss (dB) 
d 







F i gur e 



3-1 



Program "TROPQ" Module Block 



Di agr am 



36 








L = medi um-to-aperture coupling loss <dB) 
c 

L = atmospheric absorption loss (dB) 
a 

L = wavegui de/connector loss (dB) 
w 

G = transmit antenna gain (dB) 
t 

G = receive antenna gain (dB) 
r 

HG = height gain (dB) 

2. Surf ace Ret rac ti vi ty 

An adjustment to the average surface refrac— 

tivity N , refer to Figure 3-2, is made for the elevation 

o 

at each terminal site. The adjusted surface 

ref r act i vi ty , N , is [Ref. 18:p. 2-123: 



N = N exp (-0 .0 3222h ) (3-2) 

s o s 



where N = minimum monthly average refractivity 
o 

h '= average antenna height (kft) 
s 

If the surface refractivity at each site is significantly 

different, an option to calculate the respective N for 

s 

each site can be selected and an average path surface 
refractivity can be calculated. The average antenna 
height is calculated by averaging the transmit and 
receive antenna heights. 
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Figure 3-2 Monthly Sur-face Ref racti vi ty (Sea Level) 



(After Panter, Re-f . 12:p. 375) 



Effective Earth's Radius 



Changes in the surface refractivity will effect 
the radio-ray path curvature as the wave propagates over 
the earth. To accurately represent 
effective earth's radius is 

CRef. 12: p . 374D: 

a = a 1 - 0.04665 exp (0 .0055 77N ) * 

where a = actual earth's radius <6370 km) 
o 

4. Terrain Profile Plot 

. Several methods are available to plot a 

4 

troposcatter system terrain profile. The most 

fundamental method is to plot the successive path terrain 
elevations along the great circle path. Special 4/3 
earth plotting graph paper is required for this method. 
Alternatively, computer graphics techniques which obtain 
terrain information from t opogr aph i cal databases, can 
rapidly plot the profile. 

The program provides a dot-matrix printer plot. 
The terrain profile is linearly plotted by modifying the 
terrain elevations, in meters, 
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the ray-path an 
calculated as 

1 , 



to include the effect of 



the average curvature of the radio-ray path and the 

earth's surfave. Elevations, h , of the terrain are 

i 

manually obtained from topographi cal maps and tabulated 

versus distances, x , from a selected reference 

i 

terminal. The terrain data points are keyboard entered 
into the program. The modified elevations are computed 
as CRef. 12:p.380J: 




where a = effective earth's radius (Equation 3-3) 

Figure 3-3 illustrates a typical linear terrain plot. 

5. Calculation of the Angular Distance 

The terrain profile can now determine various 
path geometries. The three (3) path configurations 
considered are: 

a. Smooth Earth Horizons at Both Terminals 

b. Obstacle Horizons at Each Terminal 

c. Smooth Earth and Obstacle Horizons 

The terrain geometries are depicted in Figure(s) 3-4a, 
3-4b , and 3-4c . The respective terminal take-off angles 
are calculated for the predicted path type as 
CRef. 12:pp. 385-337 ] : 




Figure 3-3 Linear Terrain Pro-file Plot Example 
(After Panter, Ref. 12:p. 382) 
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Figure 3-4a Smooth Earth Path 
(After Ranter, Ref. 12:p. 385) 
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Figure 3-4b Near Obstacle Path 
(After Ranter, Ref. 12:p. 386) 
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Figure 3-4c Near Obstacle/Smooth Earth Path 
(After Ranter, Ref. 12:p. 385) 
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a. Type I 



Radio Horizon/Radio Horizon 



e 



V 



2 ah 



ts 



e t 



v 



2 ah 



rs 



<3-5a) 



<3— 5b ) 



er 



where h , h = transmitter, receiver terminal elevation 
ts r s 



b. Type II Obstacle Horizon/Obstacle Horizon 

_ d u 

2a (3-5 c) 



e. 



h ^+ - h K 

It ts 



et 



<9 



er 



It 



- h 

1 r rs 



lr 



lr 



2a 



( 3-5d ) 



where h ,d = transmitter obstacle el evat l on , di stance 
It It 

h ,d = receiver obstacle elevation, distance 
lr lr 

c. Type III Obstacle Hor i z on /Rad l o Horizon 



6 


= (Same 


as Equation 


3-5c ) 


et 






Q 


= (Same 


as Equation 


3-5b ) 


er 







The path type combination may be reversed to satisfy a 

Radio Hor i zon/Obstacl e Horizon configuration. Ref erring 

to Figure 3-5, the angles a and R are calculated as: 

^ oo h oo 
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Figure 3-5 Typical Path Geometry 
(After Ref. 18:p.8-6) 



46 



oo 



0 



oo 



2a 

d 

2a 



— — — + 6 <3-5e) 

, et 



h - h 

+ — S S2 * 6 



<3-5f ) 



er 



These angles are modified by correction factors, Aa 0 and 
A/3 0 to allow for the effects of a non-linear 
refractivity gradient [Ref. 12:p. 3831. The correction 
factors can be obtained from Appendix A, Figure A-7, 
however for most transhorizon, "over the horizon ", paths 
these iactars are negligible. 

The angular distance (often called scatter 
angle) is: 



6 



a 

oo 



+ 




+ A 0 O 



(3-5g) 



The ratio a oc and ft oc defines the path symmetry factor 
[Ref. 18: p . 4-71: 



S 



a oo 




<3-5h> 



The following equation will calculate the heigth of the 
intersection point of the transmit and receive antenna 
beams. This result will estimate the bottom of the 
common scatter volume (Refer to Figure 3-6) as 
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Figure 3-6 Common Scatter Volume Dimension 
(After Du Castel , Ref. 19:p. 145) 
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[Ref. 18:p. 8-83: 



S dfl 

(1 + s ) 2 



(3-5i ) 



where S = symmetry 
Du Castel , 
the dimensions of 
volume height is: 



f actor 

[Ref. 19:p 1463, approximates 
a typical common scatter volume. 



The 



d, ~ 0 , 



H = 2h 



= 0 .3h 



( 3-5 j ) 



Q , 



where 0 f = 




The maximum 



1 . 150 , 

scatter angle 
1 ong i tud i nal d 1 mensi on , 



L, 



is: 



d H 

L = = 0 . 15d 

2 h 

o 



(3-5k) 



where d = great 
The center o-f 
approximated as: 



circle path distance (m) 
gravity o-f the scatter 



volume is 



H = h + 2/3H (3-51) 

eg o 
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6 . 



Diffraction Loss 



Propagation paths having a common obstacle 
horizon, such as a mountain ridge, can be referred to as 
an obstacle gain path. It is assumed that the obstacle 
will introduce additional path attenuation. However, the 
angular distance may be reduced because of the changed 
path geometry from the obstacle. The possible reduction 
in the scatter loss may be offset by the increased loss 
due to diffraction over the obstacle. In some situations 
the common obstacle may be visible to both terminals, and 
the path loss might be less than the smooth earth path 
loss. The Internat i onal Radio Consultative Committee 
(C.C.I.R.) • has developed the following formula for 
diffraction loss relative to free-space CRef. 20:p. 1703: 




where d = transmitter to obstacle distance (m) 
a 

d = receiver to obstacle distance (m) 
b 

When the take-off angles are less than 10 degrees and 

d is greater than 2b then Equation 3-6a can be 
b a 

approximated by: 



bu 



<3-6b) 



L 



d 



20 log 



10 



2nd 




7. Worst-Hour Loss 

Seasonal annual -to-worst month path loss 
variations can be determined -from a knowledge of the 
annual changes in surface refractivity of the atmosphere 
over the path. The C.C.I.R. recommends a loss variation 
of 0.2 dB (U.S.) and 0.5 dB (Europe) per unit change of 
ref ract i ve index. 

The worst-hour median loss can be derived by 
assuming a log-normal distribution during the month. It 
was shown by Ranter, CRef. 12: p. 401 D , that on a 
log-normal distribution, the 99.9 percent point can be 



approximated by 


the value 3.1a in 


dec i bel s 


bel ow 


the 


median, where a 


is the standard 


devi at l on 


of 


the 



log-normal d i str i bu t i on . The worst-hour median loss can 
be determined as: 

Worst-Hour Median Loss = Median Annual Path Loss 
+ Difference of Annual -to-Worst-Mon th Median Loss 
+3.1 o wm 

where a = standard deviation of the worst-month 



*_j l 



wm 





Figure 3-7 



Annual and Winter Monthly Standard Deviation 
(After Ranter, Ref. 12:p. 



359 ) 



i HH 



distribution obtained -from Fiqure 3-7. 

8- Aper ture-to-Medi um Coupling Loss 

The parabolic reflector microwave antenna 
gain equation, CRef. 12:p. 103]], can be expressed in 

decibels as: 

G = 20log 10 f + 20 log -^D - 52.6 dB (3-7) 

where D = aperture diameter (feet) 

An illunimation factor of 0-54 was assumed to derive 
Equation 3-7. It would appear that Equation 3-7 depicts 
an ever increasing power gain as the antenna aperture 
area increases. However the power received by an antenna 
does not increase linearly with an increase in antenna 
diameter, D. This effect is called aper t ure-t o-medi um 
coupling loss or loss in antenna gain CRef. 12:p. 362!3. 

Aperture coupling loss in troposcatter systems 
is caused by a non-planar wavefront due to atmospheric 
l rregul ar i ti es , and a geometric effect due to the 
decrease in the scattering properties with height inside 
the scatter volume. An incoming wavefront consists of 
many plane waves, each arriving at a different angle from 
the scatter volume. If the arrival angle range variation 



is much smaller than the antenna beamwidth the wavefront 



will appear nearly plane. If the common volume is much 
wider than the receiving antenna's beamwidth a wider 
angle range will result, and the wavefront will appear 
non-pl anar . 

Basically, coupling loss can be explained as 
limited antenna pickup, as compared to the effective 
scatter volume dimensions and the antenna 3 dB 
beamwidth. Figure 3-8 compares aperture coupling loss 
results between several authors. A unique constant curve 
is presented by the C.C.I.R. CRef. 21:p. 145H. This curve 

is independent of the scatter angle and is written as: 

L c = 0 .0 7exp ^ 0 .055 (G t + G R )J (3-8) 

where G ,G = transmit, receive antenna gain (dB) 

T R 

This empirical formula gives a high coupling loss and 
will not be used in the program. The proposed empirical 
curve appears as the average of several different 
formulas and will be used as a conservative estimate for 



the aperture coupling loss. 



COUPLING LOSS 
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-o 
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Figure 3-8 Aperture-to-Medi urn Coupling Loss 
(After Levin, Ref. 21) 





9. 



Combined Free-Space/Scatter Loss 



Yeh, CRef. 221 , has derived the -following 
formula to calculate the combined free-space and scatter 
loss, in decibels: 



L - 30 log , n f + 20 log , n d + 100 + 0.2(N - 310) + 57 <3-9) 

s 10 10 s 



where f 
d 
6 



N 



operating frequency (MHz) 
great circle path distance (miles) 
scatter angle (degrees) 
surf ace refr activity 



10. Wavequi de/Connec t or Loss 

The waveguide attenuation factor was derived 
for standard rigid waveguide. At an operating frequency 
of 4.5 GHz the waveguide loss will be approx i matel y 1.25 
dB per 100 feet CRef. 18:p. 7-141. Each waveguide 

connection will introduce an additional 0.06 dB per joint 
CRef. 1 : p . 2101. 

1 1 - Absorption Loss 

Rainfall, snowfall, and fog produce atmos- 
pheric absorption loss which depends on the amount of 
moisture and on the frequency. Figure 3-9 was used to 



56 



dB per Mile 




Figure 3-9 Estimated Atmospheric Absorption 
(After Ref. 18:p. 2-43) 



/ 




estimate the rainfall attenuation on a particular 
transmission path at four (4) rainfall rates 

CRef . 18: p. 2—433 . 

C. DESIGN PARAMETERS 

1 - Carr i er-to-Noi se Ratio 

The performance of troposcatter communica- 
tions circuits are determined by the minimum acceptable 
ratio of hourly median carrier signal to thermal noise 
for a type of modulated signal. This ratio, expressed in 
decibels, is called the car r i er-t o-noi se ratio (CNR). The 
following equation is used to determine the received 
carrier power level: 



P (dBW) = P. (dBw) - L (dB) 

r t i 



(3-10) 



where P = transmitter power output (dBW) 
t 

L = total median transmission loss (Equation 3-1) 
T 

The receiver noise threshold level is written as: 



P (dBW) = -20 4 (dBW) + NF(dB) + 10 log, ^ B 
n 



(3-1 1 ) 



= thermal noise constant 
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where 



-204 dBW 



NF 



receiver noise figure (dB) 
receiver IF bandwidth (Hz) 



B 



IF 



Finally the carrier-to-noise ratio, [Ref. 10:p. 4113, 



P = receiver thermal noise level (dBW) 

n 

2- Digital Radio Link Parameters 

In digital systems the modem performance is 

usually plotted versus average bit energy, E , to the 

b 

receiver noise spectral density, N The probabilty of 

o 

error (often called bit error rate) will be determined 
from the E /N ratio. The tr ansf ormat i on of the 



1 s 



CNR ( dB) = P (dBW) - P (dBW) 



r n 



(3-12) 



where CNR 



carr i er-to-noi se ratio (dB) 



P = received power level (dBW) 
r 




E /N = CNR(dB) + 10 log 1 B - 101og 1n R (3-13) 
b o 10 w ^10 



where B 



transmission noise bandwidth (Hz) 



w 



R = transmission data rate (bit/sec) 

3. Fropability of Bit Error Calculations 

Current military multichannel communications 
use Pulse Code Modul at i on/Ti me Division Multiplexing 
(PCM/TDM) techniques to transmit digital in-formation over 
both -frequency modulated (FM) and phase modulated (PM) 
troposcatter carrier systems. The quality o-f the 

multiplexed circuits is determined by the number o-f bit 
errors that occur because o-f the channel -fading and 
multipath dispersion. The probability o-f bit errors can 
be predicted -for di-f-ferent PCM carrier modulation 
methods. Multiphase signaling and M-ary orthogonal 
signaling over a Rayleigh -fading channel are derived by 
J. G. Proakis [Ref. 14:pp. 490-499!]. 

The bit error rate (BER) -for QPSK (four-phase phase 
shift keying) and DPSK (dif f erenti al phase shift keying) is 
expressed as: 




where F = correlation coefficient 




(for coherent PSK) 
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p = 



1 + 



(for DPSK) 



where y = average received E /N per channel 
C bo 

Y. = average received E /N per bit 



b o 



h = 



L Vc 



where L = order of diversity 

j = 1 (for BPSK signaling) 
j = 2 (for DPSK signaling) 

Bit error probabi 1 i ti es are depicted in Figure 3-10 for 

two-phase and four-phase DPSK signaling with L = 1,2 and 4. 

Orthogonal signaling may be viewed as M-ary 

FSK (Frequency Shift Keying). The expression for the 

probability of symbol error (P ), derived by F'roakis, 

M 

assuming no diversity (L = 1) is: 



P M = 



M — 1 

E 

m=l 






(3-15) 



1 + m + m y 



where M = 2 (for BPSK signaling) 

M = 4 (for QPSK signaling) 

The equivalent bit error rate (BER) can be computed using: 



,k- 1 



P b = 



(3-16) 



2 k - 1 



M 
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Figure 3-10 Probability of Error (BPSK/QPSK) 
(After Proakis, Ref- 14:p- 492) 
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where k 



= log M 

The graphs of P versus E /N for M = 2, 4 and L = 1 , 

M bo 

2, 4 are shown in Figure 3-11. 



The Distortion Adaptive Receiver CDAR3 is 

currently being used in the military troposcatter digital 

radio set, AN/TRC-170. Experimental modem performance 

results are illustrated in the BER versus E /N curves, 

b o 

Figure(s) 3-12 and 3-13. The results are determined for 
three (3) different multipath delay values. 

Sunde, CRef. 24:pp. 144-2143 has derived a 

general expression for the maximum differential 
transmission delay. This equation is valid when the 
transmitting and receiving antenna beamwidths are 
di f f erent : 



5 = 



d 

2 




6 



e t 






( 3 — 1 7 a ) 



where d = path distance (meters) 

^ = scatter angle (mrad) 

6 = transmitter take-off angle (mrad) 
ct 

6 = receiver take— off angle (mrad) 

er 

The time dispersion (multipath spread) relative to the 
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Figure 3-11 Probability of Error (M-ary FSK) 



(After Proakis, Ref. 14:p. 499) 
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PERFORMANCE OF DUAL PULSE DAR 




3 8 0 12.0 18 0 20 0 24.0 28.0 02.0 



EB/NO (DB) 

(DUAL DIVERSITY L=2) 



3-12 DAR Performance - Dual Diversity 
(After Zawislan, Ref . 4:p. 4) 



re 



PERFORMANCE OF DUAL PULSE DAR 




Figure 3-13 DAR Performance - Quad Diversity 
(After Zawislan, Ref. 4:p. 4) 
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average time delay for unequal antenna beamwidths is: 



A= = ft < 3- 17b) 

2c 2(3X10 m/s) 

8 

where c = 3X10 m/sec 

In the case of equal transmitting and receiving antenna 
beamwidths, Equation 3-17b is simplified to: 

3d# 2 

\ = ( 3-1 7c ) 

16c 

For a particular order of diversity, the probability of 
error can be computed for the AN/TRC-170 troposcatter 
radiolink by evaluating a > least-squares polynomial 
equation derived from experimental data curves. 

4. Fade Margin and System Reliabilty 

The system fade margin, in decibels, is the 
difference between the "practical threshold" level and 
the median received signal level. The propagation 

reliability values for the worst fading condition, 

Rayleigh fading, can be compared with their required fade 
margins in Table II CRef l:p. 225U. The "practical 

threshold", or minimum acceptable received signal .level, 
cannot be below the FM improvement threshold 
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TABLE II 

RAYLEIGH FADING PROPAGATION RELIABILITY 




Single Hop 

Propagation Reliability 
(*/.) 

90 . 0 
99 . 0 
99 . 9 
99 . 95 
99.99 
99. 999 



Fade Margin 
(dB) 

8 

18 

28 

33 

38 

48 











[Ref. 25:p. 713. The fade margin can be evaluated by: 



Fade (dB) = P (dBm) - N (dBm) 
jo p 



(3-1B) 



where N = "practical noise threshold" 

P 

N = P (dBm) + FN improvement 
p n 

P = -174 dBm + 10 log B + NF 
n 10 IF 

W. T. Barnett and A. Vigants of Bell Tele- 

phone Labor ator i es , [Ref 25:pp. 59-603, have developed an 

empirical method to determine the nondiversity annual 

path availabilty. Barnett's procedure begins by defining 

L) as the nondiversity annual outage probability and r 
ndp 

as the fade occurrence factor: 



actual f ade _pr obab i 1 i ty 
r = Rayleigh fade probabilty 

If F is the fade margin in decibels: 

actual _f ade probabi 1 ty 
r = ~ ~ -F/10~ 

10 

For the worst month: 



r = a x 
m 




(3-19) 



where 



d 



path length (statue miles) 



f 



-Frequency (GHz) 



F = -Fade margin (dB) 

a = 4 (for smooth earth, over water, flat desert) 

a = 1 (for average terrain with some roughness) 

a = 0.25 (for mountainous terrain) 

Considering the annual fade occurrance: 



r 

yr 



br 



m 



(3—20) 



where b 
b 
b 

Finally the 



0.5 (for hot, humid coastal areas) 

0.25 (for normal , temperate or subarctic) 
0.125 (for very dry climate) 
nondiversity annual path outage is: 



U 



ndp 



r 

Y r 



10 - F/1 ° 



(3-21> 



The annual nondiversity availability percentage is: 



A = 100(1 - ) (percent) 



(3-22) 



The percentage of availability is improved 
with the use of frequency and space diversity. Figure 
3-14 is used to graphically determine the approximate 
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availability improvement for various percentages of 
frequency separation (F.S.). Figure(s) 3-15 (Dual 
Diversity) and 3-16 (Quad Diversity) provide a graphical 
method to determine the percent of path availability 
(percent of level exceeded) for different diversity 
combi ning techni ques. 

5. Pi ver si ty Requi regents 

The antenna spacing, in meters, required for 
effective space diversity has been experimentally derived 
for frequencies greater than 1 GHz, in the horizontal as 
[Ref . 11: pp . 145-146]: 

A = 0.36 (D 2 + 1600) 1/2 (meters) 

n (3— 23a) 

and in the vertical as: 

A v = 0 . 36 ( D 2 + 225) 1/2 (meters) (3-23b) 



where D = parabolic antenna diameter (m) 

A satisfactory frequency separation, in MHz, for frequency 
diversity has been dervied as: 



A f = ( 1 . 4 4 f / 



d) (D + 225 



i/: 



(MHz) 



< 3— 23c ) 
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Percent of time below level 




Figure 3-14 Approximate 
Distribution versus 
and Frequency 



Interference Fading 
Order of Diversity 
Separat i on 



(After Kef • 18: p. 2-46) 



Percent of time above level 



dB WITH REFERENCE TO The NO-OivERSITy median value 



PERCENT OF THE LEVEL EXCEEDED 




Figure 3-15 



Short— Term Fading 



(Dual Diversity 



(After Ref. 18:p. 



4 - 11 ) 



dB WITH REFERENCE TO THE NO-DIVERSITY MEDIAN VALUE 



F i gu re 



3-16 



Short-Term Fading 



(Quad Diversity) 



( Af ter Ref . 



1 8 : p . 4-13) 
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where f = transmitter -frequency (MHz) 

6 = scatter angle (mrad) 
d = path length (km) 

6. Analog Radio Link Parameters 

The analog radiolink considered will be a FM 
transmitter using -frequency division multiplexing (FDM) . 
a- Receiver IF Bandwith 

For FM modulation techniques, the re- 
ceiver IF bandwidth is computed as [Ref. 18:p. 8-161: 

BW if = 2 ( AF p + F m ) (3-24) 

where F = peak -frequency deviation (Hz) 

P 

F = maximum modulating -frequency (Hz) 
m 

The peak -frequency deviation is the product of the 
frequency modulation index and the bandwidth of the 
modulating signal. The maximum modulating frequency is 
computed as the sum of the minimum modulating frequency, 
the voice channel bandwidth, and the frequency spacing 
between multiplexed supergroups. The calculated IF 



bandwidth can now be used to compute the receiver noise 



threshold, Equation 3-11. The car r 1 er-t o-noi se ratio is 
determined by Equation 3-12. 

b. Expected Channel Noise 

According to DCA System Performance 
Spec i f i cat i ons CRef. 18:p. 3—113, the channel noise 
standard for a troposcatter link is: 



and 



L 

N(pWpO) = (16,000) 

2000 



(3-25a) 



N(dBa0) = 10log 10 (pWp) - 6 dB 



<3-25b> 



where L = path length (nautical miles) 

pUJp = picowatts psophometr i cal 1 y weight 
at, or referred to, a zero transmission level p 



The term dBa refers t 
reference noise power, 
compensate for equipment 
power that dBa is referr 
dBaO, it is required to 
this reference power 
channels, the corrected 
the expression for dBaO 



0 decibels of noise pow 
with an adjustment 

weighting. The refere 
ed to is -85 dBm. 
calculate the number o 
the signal is. For f 
reference level is - 

1 s: 



ed measured 
oi nt . 

er above a 
factor to 
need noise 
To obtain 
f dB above 
lat voice 
82 dBm and 
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dBaO = 82 - SNR 



(3-26) 



The si gnal -to-noi se ratio , SNR, must be calculated to 
compute Equation 3-26. The channel SNR may be calculated 
after the carr i er-to-noi se ratio (CNR) has been 
determined by Equation 3—11. The relationship between 
channel SNR and system CNR in a FM/FDM system is 
CRef - lsp. 2723 : 



SNR 



CNR + D. + FM. 

im lm 



P 

im 



(3-27) 



where FN 

i m 
D 

i m 
P 

i m 
L 

f 

N 



FN improvement factor (assumed 20 dB) 

diversity improvement factor (Fig» 3-17) 

preemphasis improvement factor (Fig. 3-18) 

-10 + 10 log N (dB) 

10 

number of voice channels 



where 




F 1 gu re 3- 1 7 



SNR Improvement from Diversity Technique 



(After Freeniann , Ref- l:p- 207) 



IMPROVEMENT^) 



PREEMP1IASIS IMPROVEMENT 






igure o~ 




18 Preemphasis Improvement for FDM Channe 



(After Freemann , Ref. l:p.224) 
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IV. HEIGHI^BAIN^COMPyTAIION 



A. GENERAL 

One of the purposes of this research was to develop 
a tactically practical method o-f calculating the height 
gain caused by the presence o-f an elevated troposhperic 
du.ct. The calculation must provide a result that was 
within 10 dB o-f the F’DUCT computer program prediction. 
Several related works were investigated -for a possible 
approach to the problem. Military troposcatter systems 
have introduced several constraints to the height gain 
problem. These constraints include: 

1. The transmitter and receiver terminals are both 
located outside and below the elevated duct. 

2. The operating -frequency range was between 4.5 GHz 
to 5.0 Ghz . 

3.. The geographical area of system deployment was 
limited to Western Europe (i.e. West Germany). 

A statistical analysis of elevated duct histograms 
was obtained from several radiosonde recording stations 
located throughout West Germany [Ref. 26 j. This 
information is condensed in Table III. Elevated ducts 



that had a 



percentage of occurrance greater than twenty 



percent (20 7-) over the 5 year recording period were prime 



candidates -for -further study- Eight (8) elevated ducts 
were selected as typical -for the area- Their optimum 
coupling heights ranged -from 951 to 1452 meters above the 
sur-face. The elevated duct intensities were all less 
than 6 li-units with 4 M-units being the dominant value. 

TABLE III 

ELEVATED DUCT HISTORICAL INFORMATION 



Station 
Locat i on 
(LAT/LDNG) 


Mean Optimum 
Coupling Height 
(Meters) 


Mean Duct 
In tensi ty 
(M-Units) 


Mean Du 
Thi ckne 
(Meter 


Stuttgart, FRG 
(48— 49N/09— 1 IE) 


1452 


4 


1 15 


Essen , FRG 
(51 — 23N/06— 58E ) 


1412 


3 


106 


Hannover, FRG 
(52-2BN/09-41E) 


1376 


4 


109 


Rheine, FRG 
(52-16N/ 07-25E) 


1257 


3 


1 13 


Idar — Oberstei n , 
( 49-41 N/07-19E) 


FRG 

1151 


3 


S3 


Emden , FRG 
( 53-22N/07— 1 3E ) 


1071 


4 


1 24 


Goch , FRG 
( 5 1 — 40N / 06- 1 0E ) 


1044 


6 


172 


Grei-fswald, FRG 
< 54-05N/ 1 3— 22'E ) 


951 


4 


119 



(After Ortenburger, Ref . 26, Vol . 12) 



Height gain values, in decibels, were computed by the 



PDUCT program -For each o-f the eight selected ducts. The 

-following common input parameters were entered into the 

PDUCT program -For each duct. [Ref. 71 : 

1- The transmitter site elevation was -Fixed as the 
reference at 5 meters above the surface. 

2. The receiver site elevation was increased from 
5 to 2S0 meters in 25 meter increments. 

3. The path distance was increased from 75 to 325 
kilometers in 50 meter increments. 



4. The antenna polarization was set for both hori- 
zontal and vertical polarization. 



5. The frequency of interest was increased from 

4500 MHz to 5000 MHz in 100 MHz increments. The 
frequency was changed every PDUCT program run 
for each selected duct of interest. 



6. The relative permittivity (15), conductivity 
(0.01 mho/m), and maximum mode attenuation 
(1.0 dB/km) were selected for the path character- 
i st i cs . 



B. MATHEMATICAL FORMULATION 



The PDUCT height gain results for each duct were 
investigated, and the following observations were 
discovered. The height gain values increased 
exponential ly as the receiver elevation approached the 
bottom of the elevated duct. The other significant 
observation was that for any fixed receiver elevation, 
the height gain increased linearly with increased path 



di stance- 



Two empirical methods were considered in developing 
the height gain prediction model- One method was to 
directly store the individual height gain results for 
each o-f the elevated ducts and then to interpolate a 
height gain result -from a database. This would have 
required an extensive height gain database. The chosen 
approach was to approximate height gain curves -from the 
F'DUCT results using a least-squares curve fitting 
program- A second order polynomial equation was derived 
for each frequency of interest at the initial 75 
kilometer path length for each selected duct height. 

The height gain program module consisted of eight (B) 
optimum coupling height decision regions which 
cooresponded to the selected elevated ducts- Each duct 
height region contained a height gain polynomial equation 
for each frequency of interest. For a particular 
operating frequency and optimum coupling height, a 
baseline height gain was estimated and multiplied by an 
incremental range correction factor. This range factor 
was the average differential difference in height gain 
between the successive 50 kilometer path increments. 
Figure 4—1 outlines how the height gain estimate was 



computed - 




1 .-raHmMWi wjt — bh \ 



Figure 4-1 



Height Gain Estimation Flowchart 
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c. 



RADIOSONDE DATA ANALYSIS 



Elevated ducts are described by their optimum 
coupling height, duct intensity, and duct thickness. 
This information can be determined indirectly from 
real-time radiosonde data taken along the system's path 
length. The atmospheric pressure, vapor pressure, 
temperature and relative humidity readings are used to 
calculate the refractive index, N, for various altitudes 
using the expression, CRef. 27 : p . 4 1 : 



N 

s 



77.6 
273 + T 



48.1 H P 

w 

P + 

273 + T 



(4-1 ) 



where P 
T 



P 



w 



atmospheric pressure (millibars) 
temperature (cel si us) 

saturated vapor pressure (millibars) 



H = relative humidity (7.) 

The radiosonde information can be entered into the 
program. The refractive index values are calculated for 
available altitude readings and then converted to 

modified refractive index values. At this point a 
N-profile plot can be obtained along with a radiosonde 
data listing. If an elevated duct does exist, the 
optimum coupling height, intensity, 
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and thief. ness can be 



determined directly -from the M-prof i 1 e. These duct 
descriptors must satisfy the respective limits 
established by the height gain prediction model. If the 
detected elevated duct exceeds the limits of the model, 
the height gain prediction will become inaccurate. The 
program user will be alerted of this condition. At this 
point the design engineer must decide to obtain the 
height gain from an alternative computation or neglect 



ducting effects 



V 



RESULTS 



A. GENERAL 



For validation purposes, the TROPO Program was used 
to determine the design predictions for a typical 
tropospheric scatter communications system. Table Iv 
outlines the proposed system specifications and 
assumptions. Both the terrain and radiosonde data have- 
been assumed to accommodate the program's capabilities. 
The program results are illustrated in* Figure(s) 5-1 
through 5-4, The Radiosonde Environmental Data Listing 
(Figure 5—3) identifies the occurrence of a tropospheric 
duct by detecting a "trapped" radio ray path condition. 
Other refractive bending conditions identified are 
super-ref ract .i ve (bending toward the earth's s urtacs) , 
normal (standard bending), and sub-ref ract i ve (upward 
bending). The M-Profile Plot, Figure 5-4, graphically 
verifies the program's computed duct parameters, 

B, HEIGHT GAIN COMPARISON 



Height gain prediction model results were aamDB.rad 
with the PDUCT program computed values for identical 

parameters. Height gain curves were 



elevated duct 



Height 



gain 



curves 



TABLE IV 



PROPOSED EXAMPLE SYSTEM SPECIFICATIONS 



Site: 



Tr ansmi t t er 



Recei ver 



Lat i tude: 
Long i tude : 
El evat i on : 



36 38 ' 23 "N 
06 22' 02 "W 
13 meters * 



37 09 ' 1 

05 35 ' 1 

94 meters 



Terrain Data Available: Near Obstacle Path Mode 



Radio Terminal: Military ( AN/TRC- 1 70V3) 

Diversity: Dual 

Frequency: 4560.0 MHz 

Antenna Diameter: 15 Ft. Parabolic 

Waveguide Length: 25 Ft. per Antenna 

Digital Trunk Data Rate: 2.048 Mb/s 

Transmission Bandwidth: 3.5 MHz 



Radiosonde Data Available 



Surface Duct Detected 



O K) 



plotted for two test cases: (1) a fixed optimum coupling 
height and range with a varying frequency, and (2) a 
fixed optimum coupling height and frequency with a 
varying range. Figure (s) 5~5a thru 5-5c has shown that a 
prediction error of approximately 5 dB is possible- The 
error increases as the frequency approaches the adjacent 



frequency 


i ncr ement . 


In this case height gain 


val ues 


val i dated 


at 


4700 


MHz have been averaged into 


those 


computed 


-For 


4600 


MHz to produce a shift of 


the 



prediction curve- Figure (s) 5-6 a thru 5-6c illustrate 
the results of only changing the path distance- In this 
case the greatest error detected was within 3 dB of 
PDUCT, which supports the linearity between height gain 
and range- 



TROPGSCATTER system design specifications 



SITE 

LATITUDE: 

LONGITUDE: 

ELEVATION: 



TRANSMITTER 
36.60 N 
6.37 M 
13.00 H 



RECEIVER 
37.15 N 
5.50 K 
94.00 M 



TERRAIN PROFILE TYPE: NEAR OBSTACLE PATH MODE 



TRANSMITTER TAKE-OFF ANGLE: 
RECEIVER TAKE-OFF ANGLE: 

SCATTER (ANGULAR DISTANCE): 

TRANSMITTER TAKE-OFF ANGLE: 
RECEIVER TAKE-OFF ANGLE: 

SCATTER (ANGULAR DISTANCE): 

TRANSMIT FREQUENCY: 4560.00 MHZ 



4.73 MRAD 
4.97 MRAD 
20.94 MRAD 



.27 DEGREES 
.22 DEGREES 
1.20 utsrErs 



MINIMUM RECOMMENDED FREQUENCY SEPARATION FOR QUAD DIVERSITY: 



MHZ 



AZIMUTH AT TRANSMITTER (TO RECVR) : 42. 70 (DEGREES N.) 
AZIMUTH AT RECEIVER (TO TRANS): 229. 17 (DEGREES N.) 



GREAT CIRCLE PATH: 57.91 STATUTE MILES / 93.20 KILOMETERS 

MINIMUM RECOMMENDED ANTENNA VERTICAL SEPARATION FOR SPACE DIVERSITY: IS. 5 FEET 

ESTIMATED SCATTER VOLUME BASE ALTITUDE: 417.09 METERS 



Figure 5—1 Example System Design Specifications 




Figure 5-2 Example System Terrain Profile Plot 
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MM ENVIRONMENTAL DATA LIST MM 



LEVEL 


PRESS 

(BB) 


TEBP 

(0 


VAPOR PRESS 
(MB) 


ALT 

(FEET) 


N UNITS 


N/KFT 


B UNITS 


CONDITION 


1 


1006. 0 


15.1 


15.2 


60.0 


339.7 


- 26.4 


342.6 


SUPER 


2 


1000.0 


14.2 


14.1 


281.6 


8 


15.7 


347.3 


SUB 


3 


993.0 


13.9 


15.1 


476.6 


33 d . 9 


- 11.3 


359.7 


NORBAL 


4 


962.0 


13.3 


14.3 


785.3 


sfsfyj . 4 


- 176.3 


371.1 


TRAP 


5 


972.0 


20.4 


6.0 


1071.3 


2B2.9 


26.6 


334.4 


SUB 


6 


962.0 


21.5 


S.7 


1364.9 


290.3 


.■:d o 

teU. 7 


c 5 & . s/ 


SUPER 


7 


949.0 


21.5 


6.9 


1751.3 


279.6 


.0 


363.7 


NORBAL 



TROPOSPHERIC DUCT DETECTED, TYPE: SURFACE 

OPTIMUR COUFLINB HEISHT 239 KB 

DUCT THICKNESS CCS KB 

DUCT INTENSITY (B-UNI7S) 36 



Figure 5—3 Radiosonde Environmental Data Listing 
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MODIFIED REFRACTIVITY (M UNITS) 



3SG37£t03 



.:45372E*G3 



5271tE*03 



. 3a0yj9c ,, 03 



. 76740 ^ 



. 3343i7E *03 



.34170SE*03 



. 34?043£»u3 



.3563S9E*03 



.3a3727E»03 



37,£* 



. 18IS41E-01 



. 69S239E-01 



.i2UA4£*00 



. 172904E^0 



5 .224443E+CG 



2 

O 

M 

£-• 

< 

> 

U 

u 



.2759S3E+00 



327323E+00 



. 37°063E»UQ 



.43G6G3E+00 



.48:i4:Et-;0 



3a£2£ + 00 




VE-Tic;: -IIS .'Sp CF F4GEi: "CjI-IEj FEFFACT.VITt S-uNITS’ 
rCRIISMiL ^lo; Al T 1 7u££ ."STEFS 



Figur e 



b-4 M-Prof l le Plot 



(bur f cite Dui.t ) 



.OPTIMUM COUPLING HEIGHT 
(2i9 METEKS) 



BASIC MEDIAN TRANSMISSION LOSS FACTORS 



SYSTEM LOSS FACTORS 

FREE-SPACE/SCATTER LOSS 209.4 DB 

WAVEGUIDE LOSS 7.4 DB 

CONNECTOR LOSS 4 DB 

APERTURE-TO-MEDIAN COUPLING LOSS 1.5 DE 

DIFFRACTION LOSS ilr APPLICABLE) NA DE 

F.AINFFALL ABSORPTION LOSS 3 DB 

SYSTEM GAIN FACTORS 

ANTENNA SYSTEM GAIN SB. 2 DB 

HEIGHT GAIN (IF APPLICABLE) -4.5 DB 

TOTAL SYSTEM GAIN E3.7 DB 

NET PATH LOSS 135.3 DB 

TRANSMITTER POWER 60.0 DBM 

MEDIAN RECEIVED SIGNAL -75.3 DBM 

RECEIVED NOISE THRESHOLD -91.5 DBM 

FM IMPROVEMENT THRESHOLD -SI. 5 SB 

THEORECTICAL F.F CNR 16.3 DB 

SYSTEM "ACE MARGIN 6.3 DB 



F i gur e 



i_- cr 

.J ~ vj 



System Performance Results 



SYSTEM PERFORMANCE 



SYSTEM FATH RELIABILITY 87.01 PERCENT 

EB/NO (BIT ENERGY /NOISE DENSITY) 21. S3 DB 

PROBABILITY OF BIT ERROR I353E-02 



Figure 5-5 System Performance Results (Continued) 
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Figure 5 -6a Height Gain Model Comparison 
(Frequency: 4625.0 MHz) 
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0.00 0.05 0.10 0.15 0.20 0.25 0.30 



HEIGHT GAIN MODEL COMPARISON 




H:1.352KM F:4650MHZ R:75KM 



Figure 5-6b Height 
(Frequency 



Gain Model Comparison 
4650.0 MHz) 
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Figure 5 6c Height Gain Model Comparison 
(Frequency: 4675.0 MHz) 



HEIGHT (KM) 
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v, ( **<#&*,* T^y*^**^***^ "t-Ti*?**"* h • f- 



HEIGHT GAIN MODEL COMPARISON 
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Figure 5-7a Height Gain Model Comparison 
(Range: 75 Kilometers) 
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HEIGHT GAIN MODEL COMPARISON 
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- 40.0 — 55.0 - 50.0 
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Figure 5-7b Height Gain Model Comparison 
(Range: 175 Kilometers) 
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HEIGHT GAIN MODEL COMPARISON 



HEIGHT GAIN (DB) 
H:1.257KM F:4990MHZ R:275KM 







Figure 5-7c Height Gain Model Comparison 
(Range: 275 Kilometers) 
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VI. SUMMARY 



A. CONCLUSIONS 

An interactive, computer-aided design program for 
military tropospheric scatter communications systems was 
developed. The program provides the communications 
engineer with an efficient and acccurate system 
performance algorithm that can be used in a tactical 
environment on a microcomputer. The following thesis 
objectives were accompl i shed : 

1. Both analog and digital troposcatter radio 
terminals are included as design options. 

2. The influence of elevated tropospheric ducting on 
systems installed throughout West Germany was 
considered. A height gain estimation model was 
derived using numerical results obtained from the 
PDUCT main— frame computer program. 

3. System performance predictions can be calculated 
for the newly deployed digital troposcatter radio, 
AN/TRC— 170 . 

4. Real-time radiosonde information can be analyzed 
to determine the presence of both surface and 
elevated tropospheric ducts. 

The accuracy of the height gain model was bounded by 
the range of selected duct parameters. However height 
gain results of elevated ducts which satisfied the limits 
of the model were within 5 dB of the PDUCT computer 



model- The model did not compute the height gain for 



surface-based ducts. 



B - RECOMMEND AT I ONS 

Addtional research is required and should focus 
on the following areas. 

1- The validation of the computer-derived system 
perfomance results can be conducted during 
standard operational testing. Real-time radiosonde 
data should be obtained from available tactical 
weather facilities and tropospheric ducting 
phenomena determined. Received signal strength 
data can be recorded and statistically analyzed, 

2. The height gain model should be expanded to 
included unlimited duct parameters, e.g. surface- 
based ducts and elevated surface ducts. 

3. The effects of tropospheric ducting can be studied 
further within the laboratory. Radiated energy, 
produced by 1 aser light, can be transmitted 
through a fluid medium containing 1 i ght-scatter i ng 
particles. A common scatter volume can be formed 
by mirror-like apertures at the laser sources. 

The refracti vity index of the medium could be 
controlled to simulate ducting conditions. The 
received energy could be detected and experimental 
resul ts studied- 



APPENDIX A 



The following prediction formulas were developed by 

the National Bureau of Standards CNBS) for calculating 

the long-term median basic tr ansmi sssi on loss, L , 

bsr 

[Ref. 12: p . 38911 and [Ref. 13:pp. 8-8 thru 8-1411. 

1. Long-term Median Basic Transmission Loss 



L, = 30log ln f - 20 log n D + F(0d) - F + H 
bsr ^ 10 ^ 10 o o 



where f 
d 

F< 6 d) 
F 

o 

H 



frequency (MHz) 

mean sea level arc distance (km) 
attenuation function (dB) 
scatter i ng-ef f i ci ency term (dB) 
frequency-gain function (dB) 



o 

2. Attenuation Function 

For a particular symmetry factor, S, and 

approximate surface re-fracti vi ty , N , the following 

s 

figures are used to determine the attenuation factor, 
F ( d } : 



For a surface ref ract i vi ty : 

N = 250 Refer to Figure A-i 

s 



.1 L4 



N =301 ................ Refer to Figure A~2 

s 

N = 350 ................ Refer to Figure A~3 

s 

N = 401 ................ Refer to Figure A -4 



3. Scatter i ng-Eff i cency Term 

The following equation can be used: 



F 

o 



1.086 ( 



Vs, 



(h o - 



- h 



It 



- h lr } 



( A-2) 



where 



4. 



h = height of transmi tter /recei ver antenna beam 

0 

i n ter sect i on ( km) 

h = height from obstacle elevation baseline to 

1 

the antenna beam intersection (km) 

h = transmitter obstacle elevation (m) 

It 

h = receiver obstacle elevation Cm) 
i r 

Frequency-Gai n Funct i on 

The frequency-gain function is expressed by 



as: 



H = 
o 



H (r n ) + H (r ) 
o 1 o 2 



2 



+ H (dB) 
o 



f -. M — O ) 



and 


r =41.92 


Q f h 




1 


te 




r =41.92 


6 f h 






re 



1 05 




Figure A-l Function F< 6 d) for N = 250 



(After Ref. 18:p. 8-9) 



1 06 




Figure A-2 Function F( 6 d) for N = 301 

s 

(After Ref. 18:p. 8-10) 



1 0 




Figure A-3 Function F( 0 d) -for N = 350 

s 

(After Ref. 18:p. 8-11) 



1 08 





Figure A-4 Function F( 0 d) for N = 401 

s 

(A-fter Re-f. 18:p. 8-12) 
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and 



H = 6C0.6log Yj 1 log S log q 
o 10 s 10 1.0 

where Q = scatter angle (angular distance) (mrad) 

f = frequency (MHz) 

h = transmitter elevation (km) 
te 

h = receiver elevation (km) 
re 

S = a / j 3 (symmetry factor) 

o o 

q = r / (s) (r ) 

2 1 

H (r ) and H (r ) are graphically derived using 
o 1 o 2 

Figure A-&, with 77 obtained from Figure A — 5 



1 1 0 







Figure A-5 Parameter N 

s 

Compute Ho 
(After Ref - 18: p 



Lh D Used 
o 
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Figure A-6 Frequency Gain Function 
Ref. 1 B : p . 
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Aa 0 or A 0 O n milliradians 




Figure A-7 Correction Terms -For N = 301 



(After Ref. 18:p. 6-23) 










Figure A-B Adjustment Factor for Figure A-7 
(After Ref. lB:p. 6-29) 
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It ( (HOPI.LE. 1.06) . AND. (iiOPI. GT. 1.0) ) GOTO 
IF ( (HOP'l'.LE. 1.0) .AND. (HOPT.GT. 0. 8) ) GOTO 6 
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SUBROUTINE S0RT(A,B, 

DIMENSION A ( 1) , B (1 ) 

INTEGER N,I,J 

REAL A,B,TEMEA r TEM£B 
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